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LONG-TERM  GOALS 

The  long-term  goals  of  this  work  are  to  improve  ocean  acoustic  reverberation  modeling  and  sonar 
perfonnance  predictions  in  shallow  waters  by  developing  inversion  procedures  to  estimate  seabed 
scattering  and  geoacoustic  properties  with  uncertainties,  as  well  as  investigating  the  importance  of 
various  scattering  processes.  Important  issues  include:  investigating  the  angular  and  frequency 
dependence  of  scattering  (defining  the  scattering  kernel),  determining  the  dependence  of  scattering  on 
physical  properties  of  the  seabed,  and  establishing  the  relative  importance  between  scattering  due  to 
rough  boundaries  at  the  seafloor  and  sub-bottom  interfaces  or  at  volume  inhomogeneities.  These  issues 
are  all  key  to  the  ability  to  invert  scattering  and/or  reflection  data  for  seabed  geoacoustic  and  scattering 
parameters,  and  ultimately  to  the  practical  inversion  of  active  source  reverberation  data  for  rapid 
environmental  assessment  applications. 

OBJECTIVES 

The  specific  objectives  of  this  project  are  to  carry  out  joint  trans-dimensional  (trans-D)  Bayesian 
inversion  of  scattering  and  reflection  data  to  estimate  the  in-situ  seabed  scattering  and  geoacoustic 
parameters.  This  involves  the  development  and  implementation  of  efficient  forward  modeling  and  the 
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inversion  of  several  sets  of  simulated  data  (to  ensure  the  validity  of  the  inversion  algorithm)  as  well  as 
the  inversion  of  measured  data. 
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Figure  1:  Schematic  diagram  of  the  monostatic  scattering  kernel  and  reflection 

coefficient  forward  models. 


APPROACH 

The  data  used  in  this  work  were  collected  by  Charles  Holland,  who  measured  direct-path  scattering1 
and  reflection2  data  over  a  wide  frequency  band  at  a  number  of  shallow-water  test  beds  (Malta  Plateau 
and  near  Elba  Island  in  the  Mediterranean  Sea,  and  New  Jersey  and  Scotian  Shelves  in  the  North 
Atlantic).  These  measurements  probe  the  seabed  on  an  intermediate  spatial  scale  (patch-size  radius 
-500  m  for  both  reflection  and  scattering),  which  generally  allows  a  laterally-invariant 
parameterization  of  the  seabed  and  reduces  effects  of  spatial  and  temporal  variability  in  the  water 
column  and  other  experimental  uncertainties.  While  such  data  have  been  interpreted  based  on  forward 
modeling  studies  and  supporting  geoacoustic  measurements,  rigorous  scattering  inversion  algorithms 
have  not  been  developed  previously,  nor  have  quantitative  inversion  studies  been  carried  out. 

Two  forward  models  are  used  in  this  research  to  compute  mono-static  scattering  data  and  spherical- 
wave  reflection  data,  both  of  which  are  applied  to  a  seabed  model  consisting  of  a  layered  half-space. 
Figure  1  shows  the  assumed  layered  model  structure.  The  top  (zeroth)  layer  is  seawater  and  assumed  to 
be  homogeneous  and  isotropic  with  known  properties.  Below  this  is  a  series  of  j  flat  sediment  layers, 
terminated  by  a  homogeneous  semi-infinite  basement  (j  +  1  layers  with  j  interfaces  in  all).  Layer 
properties  include  interface  depth  z  (the  lower  boundary  of  a  sediment  layer),  sound  velocity  c,  density 
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p,  and  attenuation  a.  In  addition,  the  basement  is  assumed  to  be  elastic  with  a  shear-wave  velocity  cs 
and  attenuation  as .  The  only  difference  between  the  seabed  model  for  scattering  and  reflection 
calculations  is  that  the  first  (water-sediment)  interface  is  assumed  to  be  rough  for  scattering  and  planar 
for  reflection. 

-5 

Details  of  the  derivation  of  the  scattering  forward  model  can  be  found  in  Moe  and  Jackson.  In  this 
research  the  von  Kannan  formulation4  of  the  2D  spatial  roughness  power  spectrum  (W)  is  used; 

I/k(K)  =  w2(|K|2  +  Kq)^y^2\  where  y,  wi ,  and  Kq  are  known  as  the  spectral  exponent,  spectral 
strength,  and  spectral  cutoff,  respectively,  and  K  is  the  transverse  component  of  the  incident  wave 
vector  with  magnitude  |K|  =  /cocos(0),  where  k0  is  the  wavenumber  in  the  zeroth  layer  (the  ocean). 

Reflection-coefficient  data  are  modeled  here  using  a  spherical-wave  reflection  model  instead  of  the 
more  common  and  computationally  faster  plane-wave  reflection  model,  as  the  source-receiver 
distances  in  the  experimental  procedure  are  insufficient  to  ignore  spherical- wave  effects.  Spherical- 
wave  reflection  coefficients  To  for  an  arbitrary  /V- layer  half-space  are  computed  as  in  Brekhovskikh  and 
Godin5  as  a  superposition  of  plane  waves  (the  Summerfeld  intergral). 

A  Bayesian  framework  is  applied  for  inversions  in  this  work  as  it  provides  rigorous  estimation  of 
parameter  uncertainties  as  well  as  parameter  values,  thereby  quantifying  the  information  content  of  the 
data  to  resolve  the  model  parameters.  Bayesian  inversion  is  based  on  formulating  the  posterior 
probability  density  (PPD)  of  the  model  parameters  of  interest,  which  combines  both  data  and  prior 
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infonnation.  ’  The  multi-dimensional  PPD  is  interpreted  in  terms  of  parameter  estimates  (e.g.,  the 
most-probable  model),  parameter  uncertainties  (variances,  marginal  distributions,  credibility  intervals), 
and  parameter  inter-relationships  (correlations,  joint  marginal  distributions),  which  provide  a  complete 
solution  to  the  inverse  problem.  Since  scattering  and  reflection  inversion  are  strongly  non-linear 
problems,  these  PPD  properties  are  computed  numerically  using  global  optimization  and  Markov-chain 
Monte  Carlo  (MCMC)  sampling  algorithms. 

One  of  the  goals  of  this  work  is  to  carry  out  inversions  trans-dimensionally,  a  relatively  new  approach 
in  ocean  acoustic  inversion.6  Trans-dimensional  inversion  allows  the  natural  parsimony  of  the 
Bayesian  approach  to  detennine  the  complexity  (e.g.,  number  of  seabed  layers)  of  the  estimated 
models.  It  also  allows  the  parameter  uncertainties  to  be  marginalized  over  the  number  of  layers,  hence 
accounting  for  the  uncertainty  of  the  parameterization  within  parameter  uncertainty  estimates.  Trans- 
dimensional  inversion  has  been  shown  to  model  complex  layering  structures  without  over- 

o 

parameterizing  simple  structures.  However,  trans-dimensional  inversion  is  computationally  intensive, 
and  of  particular  interest  is  implementing  population  techniques9  for  increasing  the  efficiency. 

The  inversion  of  the  two  types  of  data  (reflection  and  scattering)  will  be  conducted  both  individually 
and  jointly.  In  addition,  simulation  studies  will  be  conducted.  The  ability  to  resolve  scattering 
parameters  (spectral  strength,  spectral  exponent,  and  spectral  cutoff)  and  geoacoustic  parameters 
(sound  speed,  density,  attenuation,  and  layering  structure)  via  inversion  will  be  investigated  in  terms  of 
posterior  parameter  uncertainty  distributions,  which  quantify  the  effective  data  information  content.  In 
particular,  scattering/geoacoustic  parameter  resolution  will  be  investigated  for  inversion  of  scattering 
data  alone  versus  joint  inversion  of  scattering  and  reflection  data  (joint  inversion  is  expected  to 
overcome  inter-parameter  correlations  which  limit  parameter  resolution  in  scattering  inversion10). 
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Figure  2:  Left,  the  measured  scattering  data  (o),  and  right,  reflection  data  (o)  as  a  function  of  angle 

and  frequency  with  marginal  predicted  values. 


WORK  COMPLETED 

1 )  A  fixed-dimensional  inversion  of  the  measured  scattering  data  was  conducted;  the  results  have 
been  presented  at  the  2011  meeting  of  Acoustic  Society  of  America.11  The  inversion  assumed 
only  one  sediment  layer  above  a  fluid  basement  and  that  the  spectral  cutoff  was  fixed  at  zero. 
Geoacoustic  parameter  estimates  were  found  to  be  unresolved  (their  PPDs  were  approximately 
uniform).  An  informative  prior  distribution  constraining  the  relationship  between  sediment 
sound  speed  and  density  was  developed  based  on  a  large  collection  of  measurements  reported 
by  Hamilton.  ~  Geoacoustic  parameter  estimates  were  found  to  be  adequately  resolved  (their 
PPDs  had  a  single  clear  peak)  when  the  joint  prior  was  incorporated  in  the  inversion.  The 
scattering  parameters  were  well  resolved  in  both  cases,  with  estimated  values  consistent  with 
expectations  for  the  site  and  reasonably  small  uncertainties. 

2)  A  joint  trans-D  scattering/reflection  inversion  algorithm  has  been  developed  and  applied  to 
measured  data.  During  the  inversion  it  was  assumed  that  the  basement  was  inelastic  and  that 
data  residuals  were  independent.  Both  of  these  assumptions  were  found  to  be  incorrect. 
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Figure  3:  The  ID  (top-row)  and  2D  (bottom  row)  marginal  posterior  probability  distributions  of  the 
scattering  parameters  (scaled  so  each  distribution  has  the  same  area). 


3)  A  new  mathematical  formulation  of  the  prior  distribution  for  the  location  of  interfaces  between 
sediment  layers  was  developed.  This  formulation  allowed  for  the  explicit  evaluation  of  the  prior 
distribution  (which  is  not  possible  in  other  trans-D  inversions6’13),  which  in  turn  allowed  an 
elastic  basement  to  the  incorporated  into  the  inversion.  Ajoint  trans-D  inversion  was  carried 
out  which  assumed  errors  (residuals)  are  correlated  at  all  frequencies  and  modeled  by  first- 
order  autoregressive  AR(1)  process  with  unknown  coefficients.  The  results  were  presented  at 
the  2012  European  Conference  on  Underwater  Acoustics.  Some  highlights  are  shown  in  the 
results  section. 

4)  A  simulation  study  on  the  impact  of  data  residual  models  was  conducted.  The  study  considered 
inversion  of  simulated  data  with  significant  residual  correlation  only  at  low  frequencies. 
Inversions  were  conducted  assuming  (incorrectly)  no  residual  correlation,  correlation  at  all 
frequencies,  and  using  a  trans-D  auto-regressive  model.  A  paper  has  been  written  and  submitted 
to  JASA  and  results  will  also  be  presented  at  the  2012  Canadian  Acoustics  Association  meeting. 
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Figure  4:  Marginal  posterior  profiles  for  geoacoustic  properties  from  joint  scattering  reflection 
inversion.  Plot  boundaries  correspond  to  prior  bounds. 


RESULTS 

This  section  presents  the  joint  trans-D  scattering  and  reflection  inversion  results  from  measured  data. 
The  scattering  data  are  at  frequencies  of  600,  900,  1200,  1800,  2400  and  3600  Hz  and  an  angular  range 
of  5-24°.  The  reflection  data  have  a  similar  frequency  range  (630,  800,  1000,  1600,  2500  and  4000 
Hz)  and  an  angular  range  of  20-85°.  Figure  2  shows  the  scattering  data  (left)  and  reflection  data  (right) 
with  the  marginal  predicted  data  from  the  inversion.  There  is  reasonable  agreement  between  the 
measured  and  predicted  data;  however,  there  is  residual  error  correlation  present  for  many  of  the 
frequencies. 

The  marginal  probability  distributions  of  the  scattering  parameters  are  shown  in  Fig.  3.  These 
distributions  indicate  that  y,  wi,  and  K0  are  well  resolved  with  reasonable  estimated  values. 

Additionally  the  shapes  of  the  distributions  appear  to  be  approximately  Gaussian  and  log-Gaussian, 
respectively  (their  theoretical  distributions  if  all  other  parameters  are  known).  In  addition  the  parameter 
values  indicate  an  RMS  roughness  of  ~0.1  m  and  a  correlation  length  of  ~0.7  m  both  of  which  are 
consistent  with  the  expected  marine  environment. 
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Figure  5:  Marginal  posterior  distribution  of  basement  geoacoustic  parameters. 


The  marginal  posterior  profiles  of  the  geoacoustic  parameters  are  displayed  in  Fig.  4.  Overall,  3-5 
sediment  layers  are  resolved/supported  by  the  data.  The  seabed  sound  speed  profile  appears  better 
constrained  within  its  prior  bounds  than  the  density  and  attenuation  profiles.  The  sound  speed  and 
density  of  the  uppermost  layer  is  consistent  with  sand  (known  to  represent  the  surficial  sediments).  The 
high  sound  speed  and  density  below  5-m  depth  is  likely  evidence  of  a  limestone  basement,  which  is 
known  to  exist  in  the  region.  Figure  5  shows  the  ID  marginal  posterior  distributions  of  the  basement 
parameters.  There  is  better  resolution  of  cs  than  cb.  This  results  because  the  contrast  between  the 
sediment  sound  speed  and  the  basement  shear-wave  speed  is  smaller  than  between  the  sediment  and 
the  basement  sound  speeds. 

IMPACT/APPLICATIONS 

The  performance  of  naval  sonar  systems  in  shallow  waters  is  strongly  influenced  by  acoustic 
interaction  with  the  bottom,  and  therefore  knowledge  of  seabed  scattering  and  geoacoustic  parameters 
and  their  uncertainties  is  required  to  predict  and  optimize  sonar  performance.  Bayesian  inversion 
methods  offer  a  powerful  framework  for  parameter  extraction  and  uncertainty  estimation,  thereby 
quantifying  the  geoacoustic  information  content  of  the  data.  The  proposed  inversion  methodology  has 
been  applied  previously  to  reflection  data,  and  the  current  project  combines  this  with  monostatic 
scattering  inversion  to  constrain  seabed  scattering  parameters  and  investigate  scattering  mechanisms. 
Inversion  results  to  date  for  both  measured  and  simulated  data  indicate  reasonably  good  ability  to 
resolve  Seafloor  scattering  parameters  from  acoustic  measures. 

RELATED  PROJECTS 

1)  Quantifying  Geoacoustic  Uncertainties  and  Seabed  Variability  for  Propagation  Uncertainty 

(Award  Number.  N000140910394).  This  project  involves  Bayesian  inversion  of  reflection  data  to 
investigate  the  uncertainty  and  variability  of  seabed  geoacoustic  parameters,  and  the  effects  on 
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acoustic  propagation.  Inversion  methodologies  similar  to  that  in  the  present  project  are 
developed  and  applied. 

2)  Bayesian  Ambient  Noise  Inversion  for  Geoacoustic  Uncertainty  Estimation  (Award  Number: 
N000141110214).  This  project  involves  trans-dimensional  Bayesian  inversion  of  oceanic 
ambient  noise  for  seabed  geoacoustic  parameters.  The  ambient  noise  is  processed  to  produce 
bottom-loss  as  a  function  of  angle  and  frequency,  which  is  similar  to  the  reflection  data  inverted 
in  this  project.  Similar  trans-dimensional  inversion  procedures  are  also  applied. 

3)  Automated  geoacoustic  inversion  and  uncertainty:  Meso-scale  seabed  variability  in  shallow 
water  environments  (Award  Number:  N000140910394).  The  project  advances  and  carries  out 
geoacoustic  inversions  in  2D  and  3D  shallow-water  environments.  The  resulting  geoacoustic 
models  will  represent  benchmarks  for  meso-scale  variability  and  uncertainty  estimation,  and  also 
allow  the  study  of  compressional-  and  shear-wave  dispersion  and  attenuation-frequency 
dependence.  Bayesian  hierarchical  models  and  trans-D  inversions  allow  increasingly  automated 
data  analysis  in  challenging  shallow- water  environments. 
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